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A general method is developed for determining the inductance and characteristic impedance of 
uniform transmission lines. A non-uniform current distribution is allowed in the transverse plane. 
The system is represented by a matrix equation which can be programmed for computer solution. The 
correct inductance and impedance are obtained as the result of a simple limiting process. The method 
is applied to one particular geometry, a four-tape stripline system. Results are given for the induc- 
tance, resistance, and current distribution as functions of frequency and resistivity for a particular 
geometry. A method for extending the results to strip lines with proportional dimensions is developed. 
An accuracy of one part in 10 5 was found to be feasible for the determination of the inductance per 
unit length. 
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1. Introduction 



The determination of the inductance, capacitance, 



impedance of coaxial systems has 
of numerous reports. The system 
is the special case of a strip- 
consisting of four vanishingly thin 



or characteristic 
been the subject 
considered here 
transmission line 
parallel tapes. 

The approach used in this paper is to determine the 
current distribution in a transverse plane and therefore 
the transverse magnetic field configuration. Knowing 
the current distribution, the effective inductance per 
unit length is found directly. The characteristic 
impedance of a lossless line can then be determined 
from the inductance. 

Almost all previous authors have used the approach 
of calculating the capacitance per unit length, C, after 
determining the configuration of the transverse electric 
field. The result is then used to calculate the charac- 
teristic impedance of the lossless line, using the 
relationship 



Zo=\lvC 



(1) 



where v is the velocity of propagation in the line. 
Several applications of this approach are listed in the 
references [1— 4]. 1 The results of Cohn [2] and Bates 
[3] give analytic results which are useful design tools. 
The approach used in this paper, while nonanalytic, 
is quite general and without any geometric limitation 
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in the transverse plane. The effect of finite conductor 
losses and frequency dependence can also be included. 
It is not necessary to assume any distribution or magni- 
tude for the interdependent variables as this approach 
determines all necessary information simultaneously. 



2. Consideration of Finite Losses 

\i the line has finite losses, the characteristic im- 
pedance and propagation constant, y, can be calculated 
from 



R ±J(oL 



V G -\-jo)C 

y = a+j(3= V(R+j<oL)(G+j<oC) 



(2a) 
(2b) 



where R and G are the series resistance and shunt 
conductance per unit length of the line and a and fl are 
the attenuation and phase constant of the line. For a 
low-loss line, (2a) and (2b) can be closely approximated 
by 



Zo = vL 



i+i 



\2Q C 2qJ\ 



a z 



R GvL 

2vL 2 



P = co 



(3a) 



(3b) 



(3c) 
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where Ql — ojIJR and Q c = ojCIG. Terms involving 
Q 2 have been neglected. It is also assumed that 

o> R s 

-— > 1, and — <^ 1 

which must be satisfied for transmission line theory 
to be valid [5] where o> and e c are the conductivity and 
dielectric constant of the conductor, R s is the skin- 
effect surface resistivity of the conductor, and 77 is the 
intrinsic impedance of the dielectric. 

Both R and L are functions of frequency and resis- 
tivity as well as geometry for lines with finite loss. 
The inductance is therefore a function of resistance 
so that R must also be calculated to obtain accurate 
values of L. A computer method for calculating R 
and L as a function of frequency, resistivity, and 
geometry is developed in this paper. Values of v and 
G are readily obtained from properties of the dielectric 
material in the line. These values together with R and 
L can be used to obtain the attenuation constant from 
(3b). If the contribution of Ql and Q c in (3a) can be 
ignored, the characteristic impedance can be obtained 
from the simple relationship 



report by Hoer and Love [6] gives exact equations for 
the inductances of rectangular conductors carrying a 
uniform current. Using these equations to calculate 
the A/y, the exact d-c inductance of any rectangular 
coaxial system can be obtained. 



4. A-C Resistance and Inductance 

When an alternating voltage is applied to the sys- 
tem, the current is no longer uniformly distributed 
throughout the cross section of the conductors. The 
inductance as well as resistance then differs from 
the d-c values. 

Consider a short length of a long uniform line, such 
as shown in figure 2. Writing V2 in terms of V\ gives 



V 2 =Vt cosh yAl-hZo sinh 7A/. 



(5) 



Choose A/ such that 
approximated by 



|yA/| < |. Then (5) can be 



V 2 =V 1 -I l Z yM 



which is analogous to (1). 



3. D-C Inductance 

Figure 1 shows the cross section of a four-tape strip- 
line. The conductors marked 1 and 2 are to act as a 
single outer conductor while conductors 3 and 4 act as 
a single inner conductor in a go-and-return circuit. 
For conductors with finite resistance the current will 
be uniformly distributed throughout the outer and 
inner conductors at zero frequency. The total induc- 
tance of the system, J?, can be calculated from 



1 4 4 



(4) 



where A/y is the mutual inductance between tape i and 
j, and Mu is the self inductance of tape i. A recent 



V-x-Vi 

M 



=IZ y. 



(6) 



Substituting (1) and (2) into (5) gives 



AV 



Vx-Vt 

A/ 



= I(R+j<aL). 



(7) 



Our object now will be to calculate the total current, 
/, for some arbitrary value of AF/A/ and from this 
result obtain R and L from (7). With R and L known, 
the transmission line constants may be obtained 
from (3). 

To determine the total current in the stripline, each 
t ape is mathematically divided into a number of 
smaller, parallel sections as shown in figure 3. An 
equivalent circuit of the transmission line then looks 




v, 


H 


v 2 


-« — I 









\-At—\ 



Figure 1. Cross section of a four-tape stripline. 



FIGURE 2. Short length of a long uniform transmission line. 
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Figure 3. Cross section of the stripline showing how the 
sections are labeled. 




FIGURE 4. Equivalent circuit for the stripline with the tapes 
divided into small parallel sections. 



like that in figure 4. The width of each section is 
chosen small enough so that the current density may 
be considered uniform throughout that section. The 
appendix contains dc inductance equations which 
can be used to calculate the self, and the mutual 
inductance per unit length of these sections. The 
resistance per unit length of each section will be the 
dc resistance. 

Although the method of subdivision is almost com- 
pletely arbitrary, one system will be explained, which, 
although strange in appearance, will simplify later 
expressions considerably. Each tape is divided into 
2/i sections labeled in such a way that the index 
increases as the y axis is approached from any edge 
as shown in figure 3. The widths of the sections need 
not be equal and indeed will not be in later calculations. 
Since the geometry has symmetry about the origin and 
both axes, this symmetry is retained in subdividing the 
tapes. Symmetry produces the following relations 
between the currents in the different sections. For 
the outer tapes, 

Ij =z l2n+j = /47i+j == /6n+J9 j = 1» 2, . . . , n. 



For the inner tapes, 

l ii t j = «3n * j — *5W f j = 'In + j* J ~ 1 * 2, . . . , Jl. 

The total voltage drop per unit length along the length 
of any one section may be written, 



cm 



I=\ 



(8) 



For the purposes of this report the d-c resistance per 
unit length, a, of each section will be used in the form 
of resistivity/area. There are Sn equations which may 
be written in matrix form as 




+jai 



Mi,, ... M u 



M 8 



/' 



(9) 



\hn j 



As a result of the current symmetry conditions, only 
2n of the equations are independent. The four groups 
A = 1 to 2/i, 2/i+ 1 to 4/i, 4/i -f 1 to 6/i, and 6/i-f 1 to Sn 
are all equivalent. 

Choosing the first group, with k=l to 2/z as the 2n 
independent equations to be solved, rewrite (9) as 





+M 



l,2n \ 



E 2ft,2fl 




(10) 



tiJ = 2 Mi, 2 pn+j, 1,7=1,2, 



2/i. 



p = 



Equation (10) can now be written in the form 
V=[R+jak4Vy. 

Since Vu and h from (8) are complex, let 
Vk = e k +i/fr, h = a k —jb k . 



(11) 



(12) 
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In terms of column matrices, 

V = E+jF, I = A-jB. 



(13) 



Substituting (13) into (11) produces two real matrix 
equations: 

E = RA + a)JtB, (14) 

F = -RB + ojJtA. (15) 

Solving for A and B produces 

A = R ~ l [(l) + (coJ/R - >) 2 ] - X [E +{<»MR- l )F] (16) 

B = R- 1 [(l) + (a>J(R- i y>]- 1 [(<oJZR- 1 )E-Fl (17) 

In these equations for A and B, E and F have not yet 
been determined. To simplify their determination 
the following matrix definitions are used: 

i// = /?- 1 [(l) + (co^/?- 1 ) 2 ]- 1 

(j) = coJtR ~ l 

T = (li, . . • , l w , — lft+1, • • • f— l2n) 

6 = (0 U . • • , On, l„ + i, • • • , lnO (18) 

Since the total current in the outer conductor must 
equal the total current in the inner conductor, we have 

n 2n 

2 ak ~ S ah -° 

A- = i A- = n + 1 



1 -\rjO. This condition applied to the above results 
produces 

e\ + e n + i = 1 

The column matrices E and F can now be written in 
terms of known matrices and one scalar unknown each: 



E=d T +e 1 T T 
F=fir T 



(21) 



2 6*- 2 6,-0 

A = 1 k = )i + \ 



or in matrix notation, using (16) and the definitions 
in (18) 



tA = t*Ii[E + <I>F] = 
TB = Tiff[<l>E-F] = 0. 



(19) 
(20) 



The scalar voltage drops in the inner sections and outer 
sections obey the following conditions: 



where the superscript T indicates the transposed 
matrix. Substituting (21) into (19) and (20) and solving 
for ei and/i produces 



ej 



(ri//(/>T r )(ri//(/)^ r ) + (t\\)6 t )(tx\)t t ) 



fi = 



(Tl//T T ) 2 + (rl//0T 7, ) 2 



(Tl//(/>fl T )(Tl//T T ) - (Tl//fl T )(Tl//</)T r ) 
(Tl//T T ) 2 + (T^(/)T 7, ) 2 



(22) 



(23) 



Every term in parenthesis on the right sides of these 
equations is a calculable scalar. The result obtained 
can be substituted into (16) and (17) and the current 
in each section can be determined. With the currents 
known, the approximate inductance and resistance per 
unit length can be calculated from (7), which in terms 
of the components of the current in each section 
becomes 



1 n 

2** 



4co 



and 



Reff — 



2> + 5> 



1 n 
4 2 ak 



(24) 



5> + 2 6 



(25) 



ff=h 



fc=l, 2, . . . , n 



5. Limiting Value of L n and R n 



fn + l— fl 



e n + 1 : 



ei 



/ = n+l, n + 2, . . . In 



The correct value of inductance and resistance per 
unit length is obtained in the limit as n becomes 
infinite. A good estimate of the correct L and R can 
be obtained from several approximations calculated 
for different values of n. The approximate values of 
L and R are plotted against n as shown in figure 5 and 6. 
The asymptotes to these curves are the values as n 
Since the value of Av/A/ is arbitrary, it is set equal to becomes infinite. 



ei + e„ + 1 +;(/i+/^i) = Ai;/A/. 
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Figure 5. Approximate inductance versus number of subdivisions, Figure 6. Approximate resistance versus number of subdivision, 
showing limiting value of inductance for infinite n. showing limiting value of resistance for infinite n. 



The curve of L versus n can be approximated by an 
equation of the form 



L n = L oc + a n — y 



(26) 



where L n is the approximate inductance obtained from 
(24) for a given n, Loo is the inductance as n becomes 
infinite, and a and tf are constants. Equation (26) can 
be solved for Loo if L n is calculated for four different 
n's chosen such that 



n equal to 4, 6, and 9, (30) gave a value of Loo=0.201352 
/xH/m. For values of n equal to 9, 12, and 16, (30) 
gave a value of Loo =0.201344. For values of n equal 
12, 15, 16, and 20, (28) gave a value of Loo= 0.201343. 
The maximum difference is 5 parts in 10 5 . The latter 
value of Loo was used to obtain the curve of (L n — Loo) 
versus n shown in figure 7. The time to run these 
sets of n on a large high speed computer was 10, 35, 
and 75 seconds respectively. 



All _ tt.3 



(27) 



Then from (26) 



Loo = L M 



(L n -L n ){L n -L u ) 

\L, n Ju n ) \L, n L n ) 



(28) 



A plot of {L n — Loo) versus n for these 4 values of L n will 
be a straight line on log — log paper only if (26) is a 
valid representation of the curve L n versus n. Thus, 
four points can be used to verify (26) as well as to calcu- 
late Loo. Once it has been shown that (26) is a valid 
representation of the curve, Loo can be calculated from 
three values of n chosen such that 



Hi. 

n 2 



Hi. 



(29) 



Then (26) gives 



Loo=L w 



(Ln -Lnf 



(L n -L n )-(L n ~ Ln ) 



(30) 



o 

\ = 

o 



40 50 



The value of Loo for the curve shown in figure 5 was FIGURE 7. Difference between approximate and limiting value of 
calculated from three sets of values of n. For values of inductance versus number of subdivisions. 
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237-679 0-67— 5 



No simple approximation was found for the curve 
of R versus n. 



6. Subdividing the Tapes 

The exponent, ^, in (26) is an indication of the rate 
of convergence of L n to its limiting value, L^. Sub- 
dividing the tapes into segments of equal width gave a 
£f of approximately 1. The variation of the current 
density across the width of the inner and outer tapes 
is shown in figure 8. This curve suggests that the 
width of the segments of the inner tapes be made 
smaller near the edges where the curve is steepest. 
One such method of subdivision is to let the width, 
u>kn+j of the kn-\-j segment be given by, 



c 



Wkn - 



(n + 2-/)*' 



A = l, 3, 5, 7. 
7=1,2, . . . 



(31) 



where C is determined by setting the sum of the widths 
equal to W\2. Thus 



C 



(» + 2-»i 



W_ 

' 2 



or 
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(32) FIGURE 8. Variation of current density in 



OF TAPES, cm 

the inner and outer tapes. 



The constant, 17, was chosen such that the rate of 
convergence was near maximum. With the outer tape 
divided into segments of equal width and the inner 
tape divided into segments having widths calculated 
from (31), an tj of 3 was sufficient to yield a £f of 
approximately 3. 



7. Typical Results 

If the line has losses, the R and L will be a function 
of the actual dimensions instead of just the ratios 
t/b, Wjb and Wolb. The R and L of a stripline having 
dimensions b = W=0.5 cm and t = 0.025 cm, was cal- 
culated as a function of frequency/, resistivity p, and 
width of outer conductor W . The thickness, T, of all 
tapes was assumed to be equal and small enough to 
permit calculation of the M\j in (10) using inductance 
equations for zero thickness tapes. 2 Machine storage 
limitations prevented including the variation of current 
density with thickness. This limits the strict applica- 
tion of these results to tapes whose thickness is of the 
order of a skin depth. For the geometry considered 
here the errors introduced even at high frequencies 
should be negligible. 
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Figure 9. Variation of inductance per unit length with q. 



The variation of inductance per unit length with 
frequency and resistivity is shown in figure 9. The 
reason for using the variable q =fTk/p is that one curve 
of L versus q is sufficient for all/, 7 1 , k, andp, provided 
that T meets the requirements just discussed. The 
scale factor A: is a dimensionless proportionality con- 
stant for striplines having proportionate dimensions 
and shape. If all dimensions, except the thickness, 
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are multiplied by a scale factor, /r, the matrix M in 
(16) and (17) remains unchanged. The matrix R~ l in 
(16) and (17) is a function of k and may be written 



R l = 




(33) 



where 




w-- 



Equations (16) and (17) can now be written, 



Tk 
A= — 
P 



(D + 



o 77. 



E + (—J^)f (34) 



B = 



Tk 



(1) + 



(^jm^y 



>Tk 



nE-F 



(35) 



Since the matrices^ and Ware independent of/, T, 
k, and p, one curve of L versus the variable q=fTk/p 
will be sufficient for all/, T, /r, and p. The same result 
is true for R/Rdc versus q. For the strip line under 
consideration, with T=10~ 4 cm, p = 2XlO -6 CI cm, 
and A:=l, a q of 5 X 10 4 hertz/fl corresponds to a 
frequency of 10 3 hertz. 

A curve of R/Rdc for the stripline versus q is shown 
in figure 10. An explanation of the shape of the curve 
can be had from an examination of the current density 
curves such as shown in figures 8 and 11. The current 
density in the outer tapes remains uniform up to about 
q = S. From ^ = 5 to about # = 5000 the current dis- 
tribution is changing. Above a q of about 5000 the 
current distribution in the outer tapes no longer 
changes. The resistance of the outer tape would 
therefore remain constant up to a q of 5, increase until 
a q of 5000 is reached and then remain constant again. 
The current density in the inner tapes remains uniform 
up to a q of 50. The current then begins to crowd to 
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FIGURE 10. Variation of resistance with q. 
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FIGURE 11. Normalized current density in outer tapes 
as a function of q. 
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Figure 12. Variation of Resistance with qfor the outer 
and inner tapes. 
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Figure 13. Variation of inductance per unit length with 
width of outer tapes. 



conductors can be divided into sections whose d-c 
self and mutual inductances can be calculated. This 
technique is used in a forthcoming paper by Brooke 
and Cruz [7] to calculate the L and hence Zo of a loss- 
less, rectangular line with a rotating center conductor. 



The authors gratefully acknowledge the assistance 
and helpful criticism of C. M. Alfred, and W. E. Little. 



9. Appendix 

The mutual inductance per unit length between two 
long, thin, parallel tapes such as shown in figure 14 
can be obtained from equation 8 of reference [6] and is 



the edges of the inner tapes. The current continues 
to move to the edges as long as the value of q increases, 
never reaching a constant distribution as the current 
in the outer tape does. The curves of the change in 
resistance of the outer and the inner tapes with q would 
look like those shown in figure 12. The sum of these 
two curves gives a curve like that shown in figure 10. 
The variation of inductance per unit length with the 
width of the outer tapes is shown in figure 13 for 
different values of q, for the dimensions given above. 



8. Summary 



02 

ad 



— \n(P 2 +x 2 )-xPTdin- lX 



E—a, E+d 



;1 (* 

-I Z7 



E+d- 



0.2 



In 2/ + 



^H/m (Al) 



where the limits, which have been retained for com- 
pactness, are substituted as follows: 

S1S3 4 

[/•(*)](*) »S(-i) l+1 /to). 



The inductance and resistance per unit length of 
strip transmission line are functions of the current 
distribution in the conductors. At frequencies where 
the current distribution throughout the conductors is 
uniform, the inductance and resistance can be calcu- 
lated exactly. For frequencies where the current 
distribution is not uniform, the conductors are divided 
into sections. From these sections the approximate 
current distribution of the conductors is calculated. 
Approximate values of the inductance and resistance 
per unit length are then calculated from the approxi- 
mate current distribution. Using high speed com- 
puters, it is not difficult to divide the conductors into 
sufficiently small sections to calculate the inductance 
to an accuracy of one part in 10 5 or better. The 
accuracy on resistance is considerably less, depending 
on frequency. Other methods of subdividing the 
tapes are being considered to improve the convergence 
of R n to its limiting value, Roo. 

The variation of the inductance per unit length with 
frequency, resistivity, and width of the outer conduc- 
tors is given for a particular geometry. Tables and 
graphs of the inductance and resistance per unit length 
for other geometries are being prepared for future 
publication. 

The technique of calculating R and L as given in this 
paper for a stripline can be applied to other lines having 
more complicated cross sections provided that the 



The origin coincides with the left edge of one tape. 
If the left edge of the second tape is in any quadrant 
other than the first as shown, one or both of the values 
of E and P will be negative. The self-inductance per 
unit length of a long thin tape is 



L, = 0.2 ln-+0.2 
a 



In 2/ + | 



/>tH/m. 



(A2) 



In both of these expressions it is assumed that the 
current density is uniform throughout the conductors. 



Figure 14. Cross sectional view of two long parallel thin tapes. 
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Note that the term 



1 0. References 



0.2 



In 2/ 



2 



(A3) 



appears in both M t and /,, as the only term involving 
length. However, if M, and L t are substituted into (8), 
the terms involving length will exactly cancel. Substi- 
tuting (Al) and (A2) into (8) gives 



V,, = r k I k +JO, J] MlJi+ja) 



0.2 In 21- 



^h (A4) 



where the M// are the expressions in (Al) and (A2) 
with the (A3) term excluded. But the last term is zero 
because 



2* 



0. 
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That is, the total current in the outer conductors is 
equal and opposite to the total current in the inner 
conductors. Therefore the Mi- in (8) may he calculated 
from (Al) and (A2) with the (A3) term excluded. 
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